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ABSTRACT 

This  report  describes  the  extension  of  the  DTNS  computer  codes  to  treat  • 
non-inertial  reference  frames.  Both  laminar  and  turbulent  flow  calculations  are  , 
performed.  Additionally,  the  applicability  of  currently  available  turbulence  models'] 
in  non-inertial  reference  frames  is  investigated.  It  is  found  that  the  typical  j 

algebraic  eddy  viscosity  and  two-equation  k—t  models  used  in  inertial  reference  i 
frames  are  not  capable  of  predicting  the  effects  of  rotation.  Two-equation  k—e 
and  algebraic  Reynolds  stress  models  that  have  been  modified,  to  account  for 
rotation,  perform  somewhat  better  but  still  have  problem  arens  and  it  is  not  yet 
known  to  what  extent  they  can  be  used  with  confidence  for  complex  flow  fields. 

ADMINISTRATIVE  INFORMATION 
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This  work  was  funded  by  ONT  under  the  6.2  Quiet  Submarine  Propulsors  Program.  I'ask 
Area  RB23C22,  Program  Element  62.^23  with  internal  DTRC  Work  Unit  Number  1-1506-160-20. 
riie  development  of  the  DTNS  computer  codes  and  methodology  was  funded  previously  by  the 
Office  of  Naval  Research. 


INTRODUCTION 

There  has  been  a  tremendous  effort  in  recent  years  to  develop  Navier-Stokes  How  solvers  for 
computing  complex  How  fields.  I’hc  results  provide  numerical  pictures  of  How  fields  by  which 
complex  Huid  dynamics  phenomena  can  be  investigated  in  more  detail  than  typically  available  with 
experiments.  As  these  Navier-Stokes  How  solvers  become  more  reliable  they  can  be  used  to 
impact  the  design  process  as  demonstrated  in  Refs.  [1]  and  [2]. 

A  primary  motivation  for  this  work  is  to  investigate  the  use  of  Reynolds  averaged  Navier- 
Stokes  How  solvers  for  studying  the  effect  of  propulsors  on  the  stability  performance  of  a 
submersible  vehicle  in  turn.  Although  this  is  in  reality  an  unsteady  flow  phenomenon  it  can  be 
modelled  as  a  body  rotating  about  an  axis  and  solved  with  the  steady  state  Navier-Stokes  etjuations 
in  a  non-inertial  reference  frame.  Another  area  of  interest  is  the  flow  through  turbomachinery  and 
propulsors.  I'hese  Hows  are  particularly  complex  because  of  the  tip  vortex  formation,  rotation  of 
the  flow  field,  and  interaction  of  blade  and  hub  boundary  layers.  Current  theories  used  in  the 
design  process  are  based  on  potential  How  assumptions,  that  are  not  strictly  applicable  for  the 
complexity  of  the  How,  or  experimental  data  bases.  Simple  body  force  models  added  to  Navier- 
Stokes  flow  solvers,  Refs.  [2-4],  have  done  an  adequate  job  of  modelling  the  gross  features  of 
turbomachiniery  flow  fields  but  do  not  provide  a  detailed  description  of  the  How  and  are  of  limited 
use  because  of  the  basic  reliance  on  potential  theories.  A  full  Navier-Stokes  solution  capability  for 
such  How  fields  can  have  a  considerable  impact  on  the  design  and  understanding  of  this  complex 
How  phenomena.  Although  the  How  is  unsteady  in  the  absolute  frame  it  is  often  steady  in  the 
rotating  frame.  By  augmenting  the  David  Taylor  Navier-Stokes  (DTNS)  flow  code  to  include  non- 
inertial  reference  frames  a  tool  will  be  available  to  begin  investigating  complex  flows  with  rotation. 

The  DTNS  Reynolds-averaged  Navier-Stokes  (R.A.NS)  flow  solvers  were  developed  under 
the  Numerical  Analysis  of  Naval  Iduid  Dynamics  Accelerated  Research  Initiative  sponsored  by  the 
Office  of  Naval  Research.  The  purpose  was  the  development  of  methods  to  facilitate  the  analysis 
of  high  Reynolds  number  Hows  in  naval  geometries  in  which  viscous  effects  can  neither  be 
neglected  nor  modelled  by  approximate  formulations.  The  codes,  including  two-dimensional 
(D'rNS2D),  axisymmetric  (D  I  NSA).  and  three-dimensional  (DTNS3D)  versions,  were  designed  to 
be  relatively  easy  to  use  for  computational  fluid  ilynamics  (CFD)  analysis.  These  codes  can  be 
applied  to  a  wide  variety  of  internal  and  external  flow  fields  as  demonstrated  in  Refs.  [1.2,5-in|- 
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An  important  aspect  of  this  work  is  an  investigation  of  which  presently  available  turbiuence 
models  are  adequate  for  the  calculation  of  flows  in  non-inertial  reference  frames.  '1  he  original 
DTNS  computer  codes  contain  the  Baldwin-Lomax  [11]  and  standard  k—e  [12]  models  but,  as 
shown  by  Speziale  [13],  these  models  are  incapable  of  accurately  predicting  rotational  effects. 
Speziale  also  points  out  inconsistencies  with  the  various  modifications  that  others  have  made  to  the 
two-equation  models  but  with  the  unavailability  of  anything  better  these  will  be  tested  for  the  flow 
in  a  rotating  channel.  The  only  other  recourse  would  be  to  include  a  second-order  closure  scheme 
such  as  that  of  So  and  Peskin  [14],  Mellor  and  Yamada  [15],  or  Launder  et  al  [16]  but  it  is  doubtful 
that  these  models,  which  require  the  solution  of  six  additional  partial  differential  equations,  would 
be  generally  applicable  in  three-dimensional  flow  fields  for  which  the  DTNS  computer  codes  were 
developed.  Additionally,  these  more  complex  models  have  their  own  problems  when  computing 
rotating  flows  as  demonstrated  by  Speziale  [17].  The  scope  of  this  work  will  be  limited  to 
conventional  eddy  viscosity  and  algebraic  Reynolds  stress  models. 


NAVIER-STOKES  EQUATIONS 

The  Reynolds  averaged  Navicr-Stokes  equations  for  steady  incompressible  flow  in  a  non- 
inertial  reference  frame  can  be  written  as  (c.f.,  Batchelor  [18]) 

dx,^‘'dxjdxj  av/'^  (1) 

This  is  the  conservative  form  of  the  equations  written  in  cartesian  tensor  notation.  Here  V-i  is  the 
rotation  rate  of  the  non-inertial  reference  frame  relative  to  the  inertial  frame,  Xi  is  the  position 
vector,  and  eijk,  is  the  alternating  unit  tensor.  The  Reynolds  stresses,  Tij,  need  to  be  modelled  for 
hirbulent  flow  and  are  zero  for  laminar  flow.  The  Navier-Stokes  equations  for  a  non-inertial 
reference  frame  take  the  same  form  as  the  equations  for  an  inertial  frame  except  for  the  last  two 
terms  on  the  right  hand  side  of  Eq.  (1)  which  are  the  centrifugal  and  Coriolis  accelerations 
respectively.  It  should  be  noted  that  Eulerian  and  translational  accelerations  are  not  included  in 
the  above  equation  because  of  the  assumption  of  a  steady  state  flow  field.  I'he  local  velocit’es  in 
the  non-inertial  reference  frame,  in  Eq.  (1),  can  be  related  to  the  absolute  velocities  in  the 
inertial  frame,  Ui,  using  the  following  relation 


t/,  ==  Ui  +  e.jknjXk  (2 

The  flow  in  a  non-inertial  reference  frame  is  also  governed  by  the  continuity  equation  which 
for  an  incompressible  flow  is 


'Ihe  continuity  equation  is  frame  invariant  and  so  takes  the  same  form  for  both  inertial  and  non- 
inertial  reference  frames.  Because  the  above  equations  (1  and  3)  have  the  same  basic  form  in  both 
inertial  and  non-inertial  reference  frames  the  DTNS  codes  can  be  directly  applied  in  a  non-inertial 
reference  frame.  The  rotation  terms  are  treated  as  source  terms  and  do  not  affect  the  basic- 
equation  solver.  It  must  be  remembered  that  the  code  will  now  solve  for  the  flow  field  relative  to 
the  rotating  coordinate  system. 


Solution  Procedure 


Only  a  brief  description  of  the  solution  procedure  is  given  here  as  details  can  be  found 
elsewhere  [5,6].  The  Navier-Stokes  equations  contain  both  first  derivative  convective  terms  and 
second  derivative  viscous  terms.  The  viscous  terms  are  numerically  well-behaved  diffusion  terms 
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and  are  discretized  using  standard  central  differences.  The  upwind  differenced  Total  Variational 
Diminishing  (T^^D)  scheme  developed  by  Chakravarthy  et.  al.[19]  was  used  for  discretizing  tee 
convective  part  of  the  equations.  The  Jacobian  matrices  of  the  convective  terms  are  used  to 
generate  eigenvectors  and  eigenvalues  for  the  system  of  equations.  The  convective  terms  are  then 
forward  or  backward  differenced  based  on  the  sign  of  the  eigenvalues.  This  produces  a  third  order 
accurate  numerical  scheme  without  any  artificial  dissipation  terms  being  added  to  the  equations  for 
stability.  The  equations  are  transformed  to  a  body  fitted  coordinate  system  and  solved  using  a 
finite  volume  procedure. 

The  artificial  compressibility  technique  of  Chorin  [20]  is  used  to  add  a  time  derivative  term 
for  pressure  to  the  continuity  equation.  This  allows  the  system  of  equations  to  be  marched  in  time 
in  an  implicit  coupled  manner  using  approximate  factorization.  The  implicit  side  of  the  equations 
are  discretized  with  a  first-order  accurate  upwind  scheme  for  the  convective  terms.  This  creates  a 
diagonally  dominant  system  which  requires  the  inversion  of  block  tri-diagonal  matrices.  The 
implicit  side  of  the  equations  are  only  first  order  accurate  but  the  final  converged  solution  has  the 
high  order  of  accuracy  of  the  explicit  part  of  the  equations. 

Laminar  Flow  Calculation 

To  test  the  DTNS  code  in  a  non-inertial  reference  frame,  the  laminar  flow  between  two 
concentric  cylinders  was  computed  with  the  inner  cylinder  rotating  and  the  outer  cylinder  at  rest, 
Fig.  1. 


This  is  analogous  to  a  turbomachinery  flow  field  with  a  hub  rotating  within  an  annulus  without  the 
complexity  of  blades.  For  this  flow  an  exact  solution  can  be  obtained  for  the  tangential  velocity 
from  the  Navicr-Stokes  equations  in  polar  coordinates  (c.f.,  Schlichting  [21])  since  the  radial 
velocity  and  all  derivatives  in  the  tangential  direction  are  zero.  However,  this  is  a  good  test  case 
for  RANS  codes,  such  as  DTNS,  which  aie  written  in  cartesian  coordinates  as  none  of  the 
velocities  or  derivatives  can  be  eliminated  and  all  the  convective  and  viscous  terms  are  tested. 
Thus,  the  X  and  Y  velocity  components  would  be  different  at  each  tangential  location  but  when 
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they  arc  conibinctl  to  torm  a  tanycnlial  velocity  compoticiit  it  slrould  be  the  same  on  eacli 
tangential  plane. 

The  typical  way  to  compute  tliis  tlow  is  to  generate  a  grid  tor  the  How  field  and  obtain  the 
.solution  in  an  inertial  frame.  1  he  tangential  velocity  on  the  inner  cylinder  would  be  set  to  some 
value  and  the  tangential  velocity  on  the  outer  cyliiuier  would  be  set  to  zero.  Note  that  for  the 
D  I  NS  codes  these  tangential  \eh^cities  at  the  boundaries  would  need  to  be  convertetl  to 
appropriate  velocities  in  the  X  and  S'  directions.  To  solve  the  problem  in  a  non-inertial  frame  the 
same  grid  can  be  used  but  irow  it  is  assumed  that  the  grid  has  an  angular  velocity  If  in  the  Z 
direction  svhere  Ril  =  U g  on  the  inner  cylinder.  In  the  non-inertial  frame  the  boundary  condition 
on  the  inner  surface  is  now  no-slip  with  U  and  i’  set  to  zero.  The  outer  cylinder  is  not  moving  in 
the  absolute  frame  and  hence  the  absolute  velocities  LI  and  V  must  be  zero  there.  With  this 
information  and  liq.  (2)  we  can  get  the  folkawing  boundary  conditions  for  the  relative  velocities  on 
the  outer  cylinder 

II  =  Uy  and  u  =  — flv. 

The  How  tiekl  was  computed  with  an  outer  cylinder  having  twice  the  radius  of  the  inner 
cylinder.  I'his  How  is  iiulcpciulent  of  Reynolds  number,  R^.(  =  U (iR /r),  and  was  computed  with  a 
Rossby  number,  Ri,{=[]R /L'  g).  of  1.  The  computed  relative  velocity  in  the  tangential  direction  is 
shown  in  I  ig.  2. 


F'ig.  2.  Relative  velocity  for  the  cmicentric  cylinders. 


As  can  be  seen  it  goes  from  0  mi  the  inner  cylinder  to  -2  cm  the  canter  cylinder.  When  this  is 
converted  to  the  absolute  frame  using  liq.  (2)  excellent  agreement  with  the  exact  solution  is 
(ditaincd  as  shown  in  Fiiz.  .V 
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Fig.  3.  Absolute  velocity  for  the  concentric  cylinders. 


I'URBULENCE  MODELLING 


Tlie  calculation  of  laminar  flow  in  a  non-inertial  reference  frame  poses  no  great  difficulty 
with  an  existing  Navier-Stokes  How  solver.  In  order  to  compute  turbulent  flow  the  Reynolds  stress 
terms  in  Hep  (1)  must  be  computed  using  an  appropriate  turbulence  model.  Full  Reynolds  stress 
closure  models  have  produced  promising  results  for  certain  rotating  flow  fields  (ie.  Launder  et  al. 
[22))  but  it  is  doubtful  that  these  models  would  be  generally  applicable  in  three-dimensional  how 
fiekls  fi  I  which  the  ITI  NS  computer  codes  were  designe  1.  'I'ne  added  computational  time  oi  six 
extra  partial  dilterential  equations  for  the  stresses,  with  questionable  degrees  of  accuracy,  is  also 
undesirable.  A  practical  alternative  is  to  model  the  Reynolds  stresses  using  the  Bt'ussinesq  eddy 
viscosity  assumption 


on,  c'// , 


Where  k  is  the  tui'nilont  kinetic  energy,  is  the  Kronecker  delta,  and  u,  is  the  eddy  viscosity. 

1  he  kinetic  energy  can  be  absorbed  into  the  pressure  term  of  the  Navier-Stokes  equations  but  it  is 
necessary  to  compute  the  eddy  viscosity  using  an  appropriate  turbulence  model.  A  variety  of 
turbulence  models  have  been  proposed  for  computing  rotating  flows  within  the  eddy  viscosity 
framework,  some  of  which  are  discussed  by  Lakshminarayana  [2,5].  Speziale  [13]  has  shown 
aicoretically  that  most  of  these  models  are  not  strictly  applicable  in  rotating  flows  bi’t  since  nothing 
better  is  currently  available  some  of  these  models  have  been  evaluated  here. 


A  variety  of  models  have  been  tested  for  the  (low  in  a  rotating  channel  as  measured  by 
.1  jlmston  ct  al.  [24],  Fig.  4.  This  is  a  fully  developed  two-dimensional  channel  flow  for  which 
centrifugal  effects  are  negligible  and  the  Coriolis  forces  dominate  the  rotating  flow  field.  The  cases 
considered  are  for  a  Reynolds  number  ot  11500  with  Rossby  numbers  of  0.069  and  0.210  and  a 
Reynolds  number  of  35000  with  a  Rossby  number  of  0.068.  These  Reynolds  and  Rossby  numbers 
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Fig.  4.  Rotating  channel  problem. 


are  based  on  the  bulk  mean  vehicity  within  the  channel  and  the  channel  diameter.  The  nonrotating 
(low  is  symmetric  about  the  channel  centerline  and  can  be  computed  quite  well  with  existing 
algebraic  and  k—e  models,  I'ig. 
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Fig.  5.  Computed  profile  for  Ro=0. 

I’he  laminar  How  for  the  rotating  cases  is  also  symmetric  about  the  channel  centerline  because  of 
the  negligible  centrifugal  forces.  The  turbulent  flow  is  nonsymmetric  about  the  channel  centerline 
and  t!ie  asymmetry  in  the  problem  is  due  to  the  Reynolds  stresses.  Therefore,  even  though  the 
tlow  itself  is  relatively  simple,  it  is  necessary  to  model  the  correct  turbulent  flow  behaviour  to 
compute  the  flow  accurately.  As  stated  by  Launder  et  al.  [22]  eddy  viscosity  models  developed  for 
inertial  reference  frames  will  produce  a  symmetric  profile  for  this  problem.  Evidence  f  this  can 
be  seen  in  Fig.  6  for  the  Baldwin-Lomax  model  at  a  Rossby  number  of  0.21.  The  profile  is 


Fig,  6.  Computed  profile  with  the  Baldwin-Lomax  model  for  Ro=0.21 


symmetric  about  the  chamiel  centerline  and  does  not  exhibit  any  effects  of  rotation  on  the  flow 
field.  Because  of  this  shortcoming  of  algebraic  eddy  viscosity  models  most  of  the  reported 
attempts  at  predicting  rotating  turbulent  How  have  relied  on  the  k—e  model. 

k—(  Model 


'I'lic  modelled  form  of  the  k—e  equations  for  an  inertial  reference  frame  can  be  written  as 
f)k  ,  Oku,  c>  ,  ,  dk  Oiii 

dt  cl\',  dx,  ^  dXi 


d.x, 


dx, 


—  e 


(5) 


dt  dx. 


d 

dx, 


{u  +  u,  Id) 


de 

dxj 


-Cl 


e  On, 
k"0  dxj 


The  values  of  k  and  c  arc  then  used  to  compute  an  eddy  viscosity  using 

u,  =  C,,-^  (6) 

Typical  constants  for  the  model  are  Ci  ==  1.44,  C2  =  1-92,  and  C ^  =  0.09.  The  k—e  equations 
have  the  same  form  as  the  Navier-Stokes  equations  with  convective  and  diffusion  terms  and  hence 
can  be  solved  with  the  same  discretization  technique.  Details  of  the  upwind  differenced  scheme  as 
applied  to  the  k—e  equations  can  be  found  in  Ref.  [6].  All  of  the  models  studied  here  have  been 
previously  applied  using  the  law  of  the  wall  boundary  conditions  of  Launder  and  Spalding  [12]  or 
the  slip  velocity  of  Kreskovsky  et  al  [25],  which  is  essentially  equivalent.  However,  these 
boundary  conditions  can  produce  substantial  errors  for  complex  three-dimensional  flows  as  alluded 
to  by  Gorski  ct  al.  [26].  Therefore,  the  present  computations  bridged  the  viscous  sublayer  using 
the  near-wall  formulation  of  Gorski  [27],  'The  k—t  equations  are  invariant  under  the 
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transformation  to  a  non-inertial  reference  frame  and  hence  tlie  equations  (5)  are  the  same  whether 
or  not  the  system  is  rotating.  As  shown  by  Raj  [28]  there  is  a  rotation  term  in  the  t  equation  but 
this  vanishes  with  the  assumption  of  isotropic  turbulence.  Hence,  the  standard  k—e  model  will 
also  produce  a  symmetric  profile  about  the  channel  centerline  for  a  Rossby  number  of  0.21  as 
shown  in  Fig.  7. 


Fig.  7.  Computed  profile  with  the  k—(  model  for  Ro=0.21 


Because  the  channel  flow  is  fully  developed  at  the  experimental  measuring  station  it  is  fairly 
straightforward  to  examine  why  the  k—e  model  produces  a  symmetric  piofile.  For  this  flow  all 
derivatives  are  zero  except  d/d\'  .  The  V  velocity  component  is  also  zero.  To  produce  an 
asymmetric  How  field  the  equations  need  an  asymmetric  term.  The  term  On  /d}'  is  asymmetric 
being  positive  below  the  centerline  and  negative  above  the  centerline  but  the  derivative  of  this  with 
respect  to  y  as  it  appears  in  the  Navier-Stokes  equations  is  symmetric.  Similarly  the  Reynolds 
stress  term  changes  sign  across  the  channel  centerline  but  its  derivative  is  symmetric.  This  is 
why  the  flow  is  symmetric  without  rotation.  For  the  rotation  to  produce  an  asymmetric  flow  field 
there  must  be  an  asymmetric  term  linked  to  the  angular  velocity,  ft.  'Fhe  centrifugal  force  term  in 
the  Navier-Stokes  equations  is  not  symmetric  across  the  channel  centerline  but  this  term  is 
negligible  for  this  How.  The  Coriolis  force  term  is  symmetric  across  the  channel  centerline  and 
cannot  produce  an  asymmetric  flow  field.  Therefore  the  asymmetric  rotational  term  must  be 
contained  in  the  turbulence  model.  There  is  no  asymmetric  term  in  the  standard  k—e  model  and 
hence  it  cannot  produce  an  asymmetric  flow.  Several  models  have  been  proposed  to  account  for 
rotation  in  the  turbulence  model  most  of  which  have  been  investigated  in  the  present  work. 

.Modified  k — £  Models 

Various  models,  two  of  which  are  investigated  here,  have  included  rotation  effects  directly  in 
the  k  and  f  equations.  Wilcox  and  Chambers  [29]  modified  the  equations  by  adding 
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to  the  kinetic  energy  equation  and 


to  the  dissipation  equation.  A  similar  approach  was  taken  by  Howard  et  al.  [30]  who  did  not  modify 
the  kinetic  energy  equation  but  replaced  the  constant  Ct  in  the  dissipation  equation  by  a  turbulent 
Richardson  number  of  the  form 

+o.4(-)2n^) 

e  (TV 


Since  no  such  terms  appear  explicitly  in  the  exact  kinetic  energy  equation  these  techniques  appear 
to  be  quite  ad  hoc.  However,  it  can  be  argued  that  adding  rotation  terms  to  the  dissipation 
equation  is  an  attempt  to  model  the  anisotropic  rotation  terms  that  appear  in  the  exact  form  of  the 
equation.  The  results  obtained  with  the  two  models  are  nearly  identical  as  shown  in  Figs.  8-10. 


Fig.  8.  Computed  profiles  with  the  modified  k—i  models  for  Re=11500,  Ro=0.069. 


The  models  do  produce  asymmetry  in  the  flow  because  of  the  du/dy  term.  Because  this  term 
changes  sign  across  the  channel  it  enhances  production  on  the  unstable  side  of  the  channel,  below 
the  centerline  where  diijd}’  is  positive  (  y  =  0  ),  and  decreases  production  on  the  stable  side  of  the 
channel  where  diijdy  is  negative.  Both  models  underpredict  the  velocity  on  the  unstable  side  of 
the  channel  for  all  Rossby  numbers.  Although  the  predictions  right  next  to  the  wall  are  quite  good 
the  predictions  do  not  compute  the  steep  rise  in  velocity  evident  in  the  experiment.  This  seems  to 
indicate  that  the  near-wall  model  used  is  adequate  and  the  discrepancies  on  the  unstable  side  are 
more  a  result  of  the  rotation  effects  in  the  k—e  models.  For  the  low  Rossby  number  cases  the 
predictions  on  the  stable  side  of  the  channel  are  very  good  out  past  the  symmetry  plane.  However, 
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there  are  considerable  differences  between  the  predicted  and  experimental  values  near  the  wall  for 
the  high  Rossby  number  case  (Ro=0.21).  This  can  be  due  to  problems  with  the  modelled  rotation 
terms  at  high  rotation  rates  or  to  improperly  accounting  for  high  rotation  rates  in  the  near-wall 
region.  The  experiments  also  indicate  there  is  a  significant  laminar  region  on  the  stable  side  of  the 
channel  for  this  high  Rossby  number  case  and  use  of  k—e  turbulence  models  may  be  inadequate 
here  in  general. 

A  recent  investigation  by  Bardina  [31]  produced  a  modified  k — e  model  for  rotation  based  on 
full  simulation  data.  The  constants  C\  and  C2  in  the  dissipation  equation  were  modified  by  a  term 
of  the  form 


-.r 

e 

This  term  is  symmetric  across  the  channel  and  hence  cannot  produce  the  asymmetric  profile  of  the 
experiment. 


Modified  Models 


An  alternative  approach  to  modifying  the  k—e  equations  is  to  solve  the  standard  k—e 
equations  but  modify  the  computation  for  eddy  viscosity,  Eq.  (6).  This  method  typically  starts  with 
the  full  Reynolds  stress  equations  which  are  then  simplified  to  produce  an  algebraic  model  for  the 
Reynolds  stresses  as  demonstrated  by  Galmes  and  Lakshminarayana  [32].  This  algebraic  model  is 
then  futher  simplified  to  obtain  a  relation  for  the  eddy  viscosity  based  on  the  angular  velocity.  The 
standard  Boussinesq  eddy  viscosity  assumption  is  still  used  for  the  Reynolds  stresses  so  only  the 
eddy  viscosity  is  modified  from  the  original  equations.  This  modification  is  introduced  into  the 
eddy  viscosity  computation  by  changing  the  constant  C^j.  One  such  variation  is  that  of  Pouagare 
and  Lakshminarayana  [33]  where 

C„  =  0.09  + 

which  reduces  to  the  standard  definition  of  C n  for  no  rotation.  A  more  involved  model  was 
developed  by  Warfield  and  Lakshminarayana  [34]  which  takes  the  form 

Cn  —  —■ t(C3  —  1)(C3-^  -f  C4  —  !)/(/)]  -f  D2) 

where 


with 


Di  =  +  2^(C,  -  1)  +  (C,  -  1)= 

O.  =  [2R,(2  -  C3)p  +4(C3  -  1)(2  -  Ci)R,-R2 


R,  = 


kn 

e 


and 


R2 


dll  /dy 

n 


P  is  the  production  of  kinetic  energy  given  by 


P  = 


du-i 

dxj 


(7) 


The  constants  C3  and  C4  were  set  to  0.6  and  1.5  respectively.  Because  of  the  change  in  sign  of 
du  f()y  across  the  centerline  these  models  produce  an  asymmetric  profile  by  increasing  C,,,  and 
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hcncc  eddy  viscosity,  on  the  unstable  side  and  decreasing  it  on  the  stable  side  of  the  channel.  Care 
must  be  taken  when  using  these  models  as  C will  increase  or  decrease  depending  on  the  sign  of 
dll  /('A'.  Near  the  centerline  dll  j'dy  approaches  zero  creating  very  large,  or  small  values  for 
which  cannot  be  allowed  to  become  negative  as  this  would  produce  a  negative  eddy  viscosity  which 
is  not  realistic.  If  C^,  becomes  too  large  there  are  also  numerical  difliculties  and  hence  its  value 
was  limited  to  the  range  0<C^(;<0.4.  The  results  for  these  two  models  are  shown  in  Figs.  11-13. 


Fig.  11.  Computed  profiles  with  the  modified  C,t  models  for  Re=11500,  Ro=0.069. 


Neither  model  performs  significantly  better  than  the  other  nor  much  differently  than  the  modified 
k—f  models  tried  earlier. 


Algebraic  Reynold.s  Stress  Models 


In  an  attempt  to  get  better  results  two  algebraic  Reynolds  stress  models  were  tried.  These 
models  compute  each  of  the  individual  Reynolds  stress  terms,  r,y  ,  with  an  algebraic  equation  so 
that  the  How  is  no  longer  assumed  to  be  isotropic.  The  first  model  investigated  is  that  of  Galmes 
and  Fakshminarayana  [32]  whose  equation  for  the  stresses  takes  the  form 


-~)+  {Pu  -  -  C3) 


CiP 


(8) 


where 


Pi,  = 


duj  diii 

dXk  '  OXk 


Rij  'IpVljj  i.^ipk  '  jk  T  ^  jpk  '  ik  ) 

and  P  is  production  given  by  Eq.  (7).  This  is  a  nonlinear  set  of  algebraic  equations  for  the  stresses 
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and  can  bo  dilliciilt  to  solve  lor  a  complex  tlow  field  but  reduces  to  a  fairly  simple  form  for  the 
two-dimensional  channel  tlow  of  interest.  The  previously  investigated  model  of  Pouagarc  and 
Lakshminarayana  [33]  is  a  simplification  of  this  model  and  basically  computes  the  shear  stress  term 
ri2  from  it.  The  inclusion  of  the  full  algebraic  model  will  differ  from  the  model  of  Ref.  [33]  only 
by  the  inclusion  of  anisotropic  effects  in  the  normal  stress  terms.  The  results  for  the  full  model  are 
shown  in  Figs.  14-16. 
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Fig.  14.  Computed  profiles  using  Algebraic  Reynolds  stress  models  for  Re=11500,  Ro=0.069. 


riicy  arc  nearly  identical  to  the  results  computed  using  the  model  of  Pouagare  and 
Lakshminarayana  [33]  indicating  that  the  normal  stresses  in  the  Navier-Stokes  equations  have  little 
impact  on  the  mean  tlow  for  this  case. 


In  all  of  the  above  models,  except  that  of  Wilcox  and  Chambers  [29],  rotational  effects 
depend  on  the  Richardson  number 


R.  = 


-2n(^  -  20) 
_ ^ _ 

(fb 


Bardina  ct  al.  [35]  showed  that  for  rotating  homogeneous  shear  flow  the  turbulent  Reynolds  stresses 
do  not  scale  with  the  Richardson  number  indicating  that  these  turbulence  models  may  not  be 
generally  applicable  to  rotating  shear  flows.  One  model  that  does  not  depend  on  the  Richardson 
number  is  that  of  Speziale  [13]  which  is  an  extension  of  his  earlier  work  [36].  This  is  another 
algebraic  Reynolds  stress  model  where  the  equations  for  the  Reynolds  stresses  take  the  form 
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Fig.  15.  C  omputed  profiles  using  Algebraic  Reynolds  stress  models  for  Re=11500,  Ro=0.21. 


Fig.  16.  Computed  profiles  using  Algebraic  Reynolds  stress  models  for  Re=35000,  Ro=0.068. 
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where  Cq  =  1.68  and 
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However,  when  the  present  flow  was  computed  with  this  model  it  produced  a  profile  which  was 
symmetric  about  the  centerline  like  the  standard  k—t  model!  To  understand  this  we  need  to  look 
at  the  contribution  the  rotation  terms  make  to  the  Reynolds  stresses.  The  normal  stresses  get 
modified  by  rotation  in  much  the  same  way  as  that  of  the  Galmes  and  Lakshminarayana  [32]  model 
with  'll  increased  on  the  stable  side  and  decreased  on  the  unstable  side.  Conversely,  "22  gets 
increased  on  the  unstable  side  and  decreased  on  the  stable  side  of  the  channel.  The  problem  is 
that  the  shear  stress  dominates  this  flow  as  we  saw  earlier  when  the  inclusion  of  normal  stresses 
made  little  difference  to  the  predictions.  With  the  model  of  Speziale  [13]  the  rotation  has  only  a 
minor  effect  on  the  shear  stress  for  the  present  flow  and  cancels  out  completely  when  the  How 
becomes  fully  developed. 


To  better  include  the  effects  of  rotation  in  the  model  of  Speziale  [13]  the  current  author 
made  a  simple  modification  based  on  the  algebraic  model  of  Ref.  [32].  In  the  full  Reynolds  stress 
eciuations  the  shear  stress  is  dependent  on  the  normal  stresses  and  this  is  modelled  by  Galmes  and 
Lakshminarayana  [32]  in  the  form 


crtjiTu 


—  Tjj)  +  additional  terms 


I'his  term  was  added  to  the  definition  of  r]2  of  Eq.  (9)  using  C  =  0.25.  The  results  of  this 
modification  are  shown  in  Figs.  14-16.  As  can  be  seen  this  modification  does  produce  the  desired 
asymmetry  in  the  flow  but  the  results  are  no  better  than  tnosc  obtained  with  any  of  the  other 
models  and  in  fact  are  slightly  worse. 


CONCLUSIONS 

The  extension  of  the  DTNS  computer  codes  to  compute  flow  in  non-inertial  reference 
frames  is  a  rather  straightforward  process.  However  to  find  a  turbulence  model  that  can 
adequately  compute  the  flow  field  is  another  matter.  Standard  algebraic  eddy  viscosity  and  two- 
equation  k—€  models  cannot  properly  account  for  the  effects  of  rotation  of  the  flow  field.  A 
variety  of  turbulence  models  have  been  investigated  across  the  spectrum  from  simple  modifications 
of  the  k—e  equations  to  somewhat  complex  algebraic  relations  for  the  individual  Reynolds  stresses. 
For  the  rather  simple  flow  in  a  rotating  channel  the  results  were  reasonable  but  all  of  the  models 
underpredicted  the  velocity  on  the  unstable  side  of  the  channel  for  all  Rossby  numbers  studied. 
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The  predictions  did  not  compute  the  steep  rise  in  the  velocity  evident  in  the  experiment.  For  low 
Rossby  numbers  the  predictions  on  the  stable  side  of  the  channel  are  very  good  out  past  the 
symmetry  plane.  Flowever,  there  are  considerable  differences  between  the  predicted  and 
experimental  values  near  the  wall  for  the  high  Rossby  number  case.  Although  it  is  difficult  to 
determine  which  model  is  best,  since  all  of  the  results  are  so  similar,  the  models  of  Wilcox  and 
Chambers  [29]  and  Howard  et  al.  [30]  seem  to  slightly  outperform  the  other  more  complex  models. 
The  above  turbulence  models  should  perform  well  enough  for  complex  flows  at  relatively  low 
Rossby  numbers,  such  as  a  body  in  turn  and  some  turbomachinery  flows,  but  they  would  be  quite 
questionable  for  high  Rossby  number  cases.  A  more  thorough  comparison  of  these  models  and 
perhaps  more  complicated  ones,  needs  to  be  done  for  more  realistic  flow  fields  before  any  firm 
conclusions  can  be  made.  Such  an  effort  is  beyond  the  scope  of  this  work. 
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